function y = normalizedNoise(x,f,s,seed)

persistent init

if isempty(init)
   randn('state',seed);
   init = 0;
end

y = cos(x*2*pi + s*randn(size(x))./f);

end